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ABSTRACT 


The failure rate A, of components with exponential 
failure times, because of the manufacturing process of the 
components, behaves as a random variable. Assuming A is 
Gamma distributed with shape parameter a Anson ie parameter 
B, the reliability estimate of the components can be 
expressed as 
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R(t) = ”一 一 一 一 一 


(1 + Bt)Š 


^ 


where а апа В аге the maximum likelihood estimators of a 
and B respectively. The accuracy of this estimate is 
analyzed when the distribution of A is discretely distributed 


rather than gamma distributed. 
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т. ZENDROPUCTION 


Let T be the life time of a component that is 
exponentially distributed with parameter A. Because of 
the manufacturing process and different logistical forces 
acting on the components, A behaves as a random variable 
шаға certain distribution G. The distribution of the 
time to fallure,under the above assumptions, is called a 
mixed exponential distribution. 

Let N components be put on test until a test time To 
is reached. Suppose that as a result of testing, 

E 


t m 


1? gr.” 


and μμ. N 


паа not failed. if Gis a Gamma distribution with share 


K were life times of components that had failed 


.„t,, were life times of components that 
parameter a and scale parameter B, then a and В сап be 
estimated by using the method of maximum likelihood. 

шиг в.в. οἱ chezreizabllity function of the components 


can be expressed as 


al 


R(t) =з з ex 
ρε 


where a is the M.L.E. of a 
Ве М. Б.Б. of B 
Computer simulations were conducted to determine the 
Sensitivity of this estimation procedure when the exponential 


distributions are mixed differently from the assumed Gamma 


aaovuribution mixture. 





ЦЕП ШЕ T TP DIII OF THE MODEL 


The model and the derivations in this section are 
discussed in detail in a paper by Myhre and Saunders [5]. 

Let T the life time of a component be exponentially 
distributed with parameter A, where A is a random variable 
with Gamma distribution G with shape parameter a and scale 
parameter В. 

The density function of A is 


о-1 -А/8 
оо asy s Ae. —_ (1) 
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The reliability of the component is 


R(t) P(T>t) 
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пней the density function of T is given by 


_ aR(t) 
dt 


f(t) 


ag 
= --------- (3) 
(Аве) т 


suppose that of N items that were put on test it was 


Observed that t.,t.,...,t, were the life times of k items 


jc k 


that had failed before 16 and let t = t 


k+1 Kk+2， o” 
were the life times of items that had not failed. The 
maximum likelihood estimators of the parameters are derived 
as follows: 


The likelihood function is 
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12012) JN E CER) (4) 
1 5 j=k+1 J 
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J op) = 
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Substituting equations (2) and (3) into (4), 
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Taking the logarithm of L(a,8), 


гэ 
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По 50122) 


k 1n(a8) - a 


Пес, 
ных 


11(1%8%,) - 1n(1+8t,) 


1-1 = 


(6) 


The maximum likelihood estimators of a and ß can be 


obtained by differentiating La 
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wo 2E. E In(1+8t,) (7) 
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1-1 
Е. (8) 
ов ізі 188% 1-1 1986, 


Suppose first that a is known. Then B, the estimate 


of 8 can be computed by setting equation (8) to zero 
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ка r 217 E (9) 





Rearranging equation (9), 
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Finally, 
Ч ЭГ”: 
+ E » - 
НЭЭ one 0 (10) 
The solution of equation (10) is B, an estimate of В. 
Note that since 
N k 
1 1 
AB) = a EZ + È - Na 
Еа 0 
is a decreasing function of ß with A(0) = k and 
lim A(ß) = -Na, there exists exactly one solution for ß 


Bo 
Ha nevera > O and N > k > 1. 


If B is known and B is positive then a could be computed 


directly by setting equation (7) to zero, 


(11) 
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If a and 8 were both unknown, a and ß could be obtained 


by solving the equations simuitaneously 
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Eliminating a from equations (12) and (13) results in an 


equation that contains only 8, say B(8) = O where 








1 N k 1 М Bt. 
Ε ο το = iIn(ltet,) X = 7 (14) 
K i=] er et 
шэг ls the solution of ВС2) = O and a can be solved 


using equation (11). 


We now proceed to examine the behavior of the function 


B(B). In particular B exists only under certain conditions. 


Note that B(0) = 0 and lim B(ß) = -N and B(ß) is a continuous 
->00 
function. Hence, B(8) = O would have a positive solution 


if its derivative at B = O is positive. Since B'(0) = O, 


1 
the sign of lim — may be examined. A condition that 
B>0 ' 
B(8) = 0 has a positive solution is that lim B'(8) > 0 
B0 


(see Figure 1). 
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Figure 1 
Graph of B(ß) 
Rewriting equation (14), 
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since it is required that lim 
^ B0 
Solution for B, it can be stated from equation (17) that 


0 со get a positive 


a sufficient condition for the existence of a positive 6 


Ns that 
k 
t ο 9 (18) 


The reliability estimate R(t) could be computed using 
equation (2) whenever a positive solution for 8 exists. 
Шаға positive solution for 8 does not exist however, R(t) 
can not be computed using this equation. In this latter 
case, the reliability could be estimated in the following 
manner. 

Assuming A is Gamma G(a,ß) distributed, then E(A) = aß. 
Setting aß = y anà substituting a = Y/B into equation (6), 


the log likelihood function becomes 
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The M.L.E. of y and B could be obtained by setting equations 
ёЕ PI) to zero, and Solving simultaneously for Y and 


8. The resulting estimator for y is 


A O 
"EM кё (22) 
Y In(1+ßt,) 
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Substituting y from equation (22) into equation (21) and 


rearranging, 








which is exactly the same as equation (14). Thus the M.L.E. 
of ß is the same as before. However whenever a positive 
solution for 8 does not exist the reliability estimate 
could be computed by treating the life times of the com- 


ponents as exponential random variables with a constant 
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failure rate y, that could be computed by taking the 


limit of equation (22) to get 


Eo ovs (23) 
EN 0) C. 
i-1 i-1 


In this case, the reliability estimate is given by 
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ОО ος (24) 
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БИН EUTER STIMULATION 


To determine the robustness of this estimation procedure, 
computer simulations were conducted for a discrete probability 
@derstribution on A. Let Ny Agar ee od be the chosen values 
of failure rates of the components with associated proba- 
bilities Ῥη»Βϱ»:»»»Βῃῃ» A total of N exponential failure 
times T's were generated using a standard exponential 
random number generator. The T's were then truncated at 
a planned test time To: If this set of T's gave any posi- 
tive solution for 8, then 8 could be computed by setting 
equation (14) to zero and solving for 8 and 2 could be 


Ebrsugmed $pomv equation (11). The reliability estimate was 


computed by using 
R(t) = —— (25) 
(1*8t)? 


If no positive 8 exists, then the exponential procedure 


^ 


was used. Y the estimate of E(A) was computed by using 


к (26) 
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where 


N: the number of components on test 


T! the planned test time 


0: 

the number of components that failed before To: 
The reliability estimate then was computed using equation 
(24). 


Since the true reliability of the component is 
m 
Е Dae (27) 


where m is the number of chosen failure rates, a measure 
of robustness of this estimate is R(t) - R(t). 
A delineation of the simulation steps is as follows: 


(1) Generate N times to failures T T5,...,T4 Where 


de N 
Ti is exponentially distributed with parameter 
A; and each T is randomly generated in accordance 


ucchesstcasted probability distribution on A. 


(2) Truncate Ti ,T5,... Ty at a test time To: 


(5 cheek if this set of Ti ,T5,.. 


B by taking a small value Во then check the value 


«Ly give any positive 


of Β(β0) from equation (1!). 

1 38 B(89) is positive, continue to step (4). 
ἘΠ B(89) = 0, then 8 - By» and continue to step (5). 
ТЕ B(89) < 0, use equation (26) and (24) to compute 


the reliability estimate then continue to step (7). 


L 








(4) Compute β from B(ß) = 0 using Newton-Raphson method. 
(5) Compute 5 ШЕШІ equacion (11). 
(6) Compute the reliability estimate R(t) using 
equation (25). 

(7) Replicate steps (1) through (6) Τ᾽ 

(8) Compute the mean value and standard deviation of 
all values of R(t), and let R(t) be the mean value 
of R(t). 

(9) Compute the true reliability R(t) by using equation 
ШЕ 

(10) The measure of robustness is R(t) - R(t). 

A second simulation was done as a comparison to the 
first one. The procedure was similar to the first one 
except that instead of truncating the failure times 
Т,,Т,,....Тұ at a test time To (see step 2), they were 
truncated at the time of the 220 fallure. In other words 
N components were put on test, and the test was terminated 
after k components had failed. 

Using this procedure however, apositive solution for β 
КЕЕ псует ехіес fork < 2. Only if k > 3 under certain 
Esnditions will a positive 8 exist. This can be shown as 
follows. 


Recall that a sufficient condition for the existence 


of a positive 8(see inequality (18)) is, 
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N > 2 N k 
за шаа! επι 
For k-1 then ty-t, ,7...-t4-t, and the inequality (18) 
becomes 
2 2 
Nt, > eNt, 


мЕНС his, of course is invalid. 


For k=2, t,,=t and the inequality (18) 


MERI 
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becomes 


2 
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t 1 


+ (N-1)t.° TIU. t 


2 
т (N-1)t, > 7 ito 


and this is also invalid. 


For k=3, t,-t the inequality (18) becomes 


Vtt 


2 2 
中 (N-2)t, > 3 (6, *t5* (N-2)t4) (t, ttt.) 


2 2 
t (N-2)t, ] > 3L2t, t5 * (N-1)t t, (N-1)t,t. ] 


and this inequality could be valid for some sets of 6, '5. 
Thus, using this’ later simulation procedure β par Wale 
possibly exist only for k E 3 
Using these two procedures, simulations were done and 


the results appear in Section IV. 


19 








Пү RESUÉTS OM THE SIMULATION 


Graphs which depict the results of the simulation 
appear on the following pages. Symbols used on the graphs 
аге: 

R(t): The reliability at time t 


R(t): асе с расе о ЕКЕ) 


N: Number of components on test 

T : The pianned test time 

m Number of replications 

А: Parameters of generated life time T's of the 
components 

р: Discrėte probability distribution for à 

k: Number of components that failed before То the 


planned test time 


PB: Number of replications that gave positive 
Ssouuvions for В 


NB: Number of replications that did not give a 
positive solution for ß 


Шеге Total number Of replications that had k failures 


To illustrate the graphs, let us observe Figure 2. This 
graph shows that five exponential life times (N=5) were 
generated via a set of à: .010, .010 with probability 
ШОО ГУ Pp: -.50, .50. By using a planned test time ia = 10, 
the reliability estimate was computed. With 50 replications 
(r=50) the average reliability R(t) was computed and was 
Shown on the graph. Further, it can be observed that for 


k = 0 (no failure occurred), NB = 28 (28 replications 


20 





occurred without any failure). It happens that 3 replica- 
tions occurred with 2 failures (Total = 3, k = 2) and all 
of these three replications gave no positive 8 (NB=3) and 
Of course no replication gave positive solutions for 8 
(PB=0). 

Another example is Figure 25 where 50 replications 
(r=50) of 10 exponential life times (N=10) were generated 
аға set of À: .005, .020 with probability density 
Ш 50, .50. It is observed that by using a planned test 
time To = 20 unit times, 22 replications (Total=22) occurred 
with 2 failures (k=2) of which 5 replications gave positive 
Solutions for 8 (РВ-5) апа the rest did not give any 
positive solution for 8 (NB=17). 

Figure 44 shows that 50 replications (r=50) of 
10 exponential life times (N=10) were generated via a 
EESTI A: .005, .010, .015, .020 with probability density 
р: .40, .20, .35, .05. The N items were tested until a 
fixed number of failures (k=5) occurred. Eight of those 
replications gave positive solutions for 8 (PB=8) and the 


rest (NB=42) gave no positive solution for 8. 
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V. CONCLUSION 


The results of the simulations show that 
l. In all cases, under the first procedure of 
simulation a positive solution for B was never obtained 


for relatively small planned test times, e.g., To = 10. 


^ 


Some positive solutions for B were obtained when To = 20. 

2. By using the second procedure (fixed number of 
failures) a positive solution for 8 never existed for 
k = 1 and k = 2. Oe Cue severe positive solutions 
Tor β obtained. The larger the value of k, the more 
accurate the computed reliability estimate will be. 

3. The result also showed that most of the reliability 
estimates have quite small error relative to the true value 
of the reliability. It was less than 1% on the average. 
The smallest relative error of the reliability estimates 
was .02% while the largest one was about 3.15%. Thus the 
proposed estimation procedure is quite accurate when the 
underlying distribution on A is a discrete distribution 
similar to those assumed rather than a continuous gamma 
distribution. In this regard the estimation procedure 


indicates some property of robustness. 
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